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Abstract 



The computation of the 0(a s ) QCD corrections to the cross sections for W + W~ r y 
and ZZ^f production in hadronic collisions is presented. We consider the case of a 
real photon in the final state, but include full leptonic decays of the W and Z bosons. 
Numerical results for the LHC and the Tevatron are obtained through a fully flexible 
parton level Monte Carlo based on the structure of the VBFNLO program, allowing an 
easy implementation of arbitrary cuts and distributions. We show the dependence on 
scale variations of the integrated cross sections and provide evidence that NLO QCD 
corrections strongly modify the LO predictions for observables at the LHC both in 
magnitude and in shape. 



1 Introduction 



The experimental precision that will be reached in cross section measurements at the 
CERN Large Hadron Collider (LHC) demands an effort of the theoretical community in 
providing accurate phenomenological predictions. NLO QCD corrections for cross sections 
and distributions have thus become mandatory, and many relevant processes at hadron 
colliders are now known to this accuracy. 

Events with multiple gauge bosons in the final state provide an irreducible background to 
many new physics searches (see, for instance, [1] for a discussion of the relevant backgrounds 
in the search of New Physics at the LHC). In addition, the triple gauge couplings involved 
in the contributing diagrams allow for restrictive tests of the gauge sector of the Standard 
Model. The process pp — > W + W~ r y + X is particularly important since it is also sensitive 
to four gauge boson couplings, namely the WWZj and WW^yy vertices [2]. 

In this paper, we present the 0(a s ) corrections for the processes 

" W+W-j" pp, pp -> viJfPi^-y + X 

"ZZ-f pp, pp - /+/r/ 2 + / 2 "7 + X. (1.1) 

Similar to previous work on triple weak boson production [3-6], we find that the QCD 
corrections are sizeable and also modify the shape of the differential distributions for many 
observables: this proves that a simple rescaling of the LO results is not adequate and a 
full NLO Monte Carlo is needed for a quantitative determination of quartic couplings at 
the LHC. We have implemented our calculation within the VBFNLO framework [7], a fully- 
flexible parton level Monte Carlo program which allows the definition of arbitrary acceptance 
cuts and distributions. 

The paper is organized as follows: in Section [2] we provide an example of the relevant 
Feynman diagrams at tree level, a short account of the strategies used to compute the 
real and virtual corrections and the various checks performed both internally and against 
other available codes. In Section [3] we show numerical results, including the scale variations 
of the LO and NLO integrated cross sections and some selected differential distributions. 
Conclusions are given in Section HI 

2 The calculation 

We consider the processes (11 .ip up to order a s a 6 in the limit where all fermions are 
massless. Among the LO diagrams (110 in the u W + W~ r j n case, 336 in the u ZZj" case) 
we can distinguish three different topologies: they correspond to the cases when 1, 2, or 3 
vector bosons are attached to the quark line (respectively case a), b), or c) of Figfj], where 
examples for u W + W~ , y" production are shown. The single and double vertex topologies 
of FigJT] a) and b) also exist for "ZZ7" production since we include off-shell effects in our 
calculation, i.e. the photon or a Z-boson can be radiated off a final state lepton.) 
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Figure 1: Examples of the three topologies of Feynman diagrams contributing to the process 
pp ^W + W~'j+ X at tree-level. 



In all cases, the leptonic decays of the W and Z bosons and combinations of subgraphs 
corresponding to decays of the same virtual W, Z, or photon are factorized in the form of 
leptonic tensors and computed independently from the rest of the cross section, in analogy 
with the procedure in [8] : this reduces computational time since the same decay current may 
appear in many Feynman diagrams and in different subprocesses and can, thus, be computed 
once for each phase-space point and saved for later use. For the computation of the matrix 
elements, we use the helicity method introduced in [9]. 

The cross section at order a s contains contributions from real emission and virtual cor- 
rections from the interference of one-loop diagrams with the Born amplitude. We use the 
Catani-Seymour dipole subtraction method [10] to handle the cancellation of infrared diver- 
gences between real emission and virtual corrections. Factorization of initial state collinear 
singularities into the parton distribution functions leads to additional "finite collinear terms" . 

The NLO real corrections are generated by the emission of a gluon off the quark line 
or by the emission of a quark through a g — > qq splitting in the initial state. The strategy 
of computing leptonic tensors separately proves to be particularly useful in this case, since 
we have a total of 272 (896) different Feynman diagrams contributing to the L W + W~ , y" 
( ) process when attaching an additional gluon to the Born diagrams. 

In the case of real emissions, the definition of an isolated on-shell photon requires further 
discussion. Photon emission collinear to a massless quark results in an additional singularity, 
which needs to be regularized. The simple rejection of events containing partons inside a 
cone of fixed width drawn around the photon is not a permissible choice since it will suppress 
the phase space for soft gluon emissions and spoil the cancellation of infrared divergences. 
In this paper we avoid the complications of quark fragmentation to photons by using an 
effective method which has been proposed by Frixione [11]. The precise definition of this 
photon isolation cut is discussed in Section El 

The NLO virtual corrections come from the interference of leading order diagrams with 
one-loop diagrams obtained by attaching a gluon to the quark-antiquark line. In the case 
where only one vector boson is attached to the quark line, i.e. topology a), one only has 
vertex corrections, while diagrams belonging to topologies b) and c) give rise to boxes and 
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pentagons. The tensor coefficients of the loop integrals have been computed a la Passarino- 
Veltman [12] up to box level, while the Denner-Dittmaier reduction procedure [13] has been 
used for the pentagons, which avoids the problem of small Gram determinants occurring in 
pentagons with planar configurations of the external momenta. Denner-Dittmaier reduction 
is now fully implemented in the public version of the VBFNLO code. The overall contribution 
from virtual corrections can be written as 



M V = M V + ^-C F r(l + e) 



2 3 4 



M B , (2.1) 



where Mb is the Born amplitude and Q is the partonic^ center-of-mass energy, i.e. the 
invariant mass of the final state VVj system. The term My is the finite part of the virtual 
corrections to 2 and 3 weak boson amplitudes as in FigJTb) and c) which we call boxline and 
pentline contributions in the following, so named since boxes and pentagons constitute the 
most complex loop diagrams, respectively. 

We have verified analytically and numerically that in order to obtain the above factor- 
ization formula of the infrared divergences against the born amplitude, the transversality 
property of the photon (e(p 7 ) • p 7 — 0) must be used, otherwise, additional IR singularities 
which are not proportional to the born amplitude do appear. After checking that there are 
no additional 1/e terms due to the massless on-shell photon, the analytical structure of My 
can be obtained from the massive weak boson case considered in [5]. Nevertheless, we have 
recomputed equivalent pentline and boxline contributions for massless particles and found 
that they agree at the double precision accuracy level with those of [5], once the proper 
scalar integrals are used. 

A further reduction of computing time is obtained for u W + W~ , y" production by imple- 
menting a trick already used in [4,5,8]. We write the polarization vector of the vector bosons 
as 

e v = x v q v + i v (2.2) 

where qy is the momentum of the massive boson (V = W^). A pentagon contracted with 
one of the external momenta qy can always be reduced to a difference of box diagrams. 
Thus we effectively isolate a true pentagon contribution and shift the rest of the 5-point 
contribution to the less time-consuming boxes. The remainder iy is chosen in such a way 
that 

ev • + qv 2 ) = , (2.3) 

i.e., the time-component of the shifted polarization vector is zero in the center-of-mass system 
of the V-pair. This choice yields particularly small pentline contributions which, hence, 
can be determined with lower Monte Carlo statistics. For "ZZ7" production, the pentline 
contributions group into gauge invariant subsets which are invariant under the replacement 
of Eq. (12.21) . Thus the trick does not speed up the code but provides a consistency check for 
the calculation. 

We have performed a number of checks on our final results. First, we compared all tree 
level matrix elements (including the ones for real emission corrections) against amplitudes 
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generated by MadGraph [14]: the agreement is at the level of machine precision (10~ 15 ). 
Second, we have compared the LO cross sections both for pp W + W~ r y + X and pp — > 
ZZ^ + X against HELAC [15] and found agreement at the per mill level. Third, we have 
compared integrated tree level cross sections for pp{pp) — > W + W~~^(j) + X and pp{pp) — > 
ZZ'j(j) +X against MadEvent and Sherpa [16] finding agreetment at the per mill level. Fi- 
nally, we also performed gauge tests based on Ward identities for the NLO virtual corrections 
at different levels of complexity. 



3 Results 

Numerical results have been obtained through the implementation of our calculation into 
a NLO Monte Carlo program based on the structure of the VBFNLO code [7]. This and other 
processes with a real photon in the final state will be included in the public version of VBFNLO 
in a future release. 

The value of the total cross section for "ZZ7" production, including leptonic Z-decays, 
turns out to be less than a femtobarn at LHC energies: since it will probably be too hard to 
detect, at least in an early luminosity run, in this section we focus on the pp — > W + W~ r y+X 
process. 

For the electroweak parameters we use the W and Z boson masses and the Fermi constant 
as input. From these we derive the electromagnetic coupling and the weak mixing angle via 
tree level relations, i.e. we use 

m w = 80.398 GeV sin 2 (9 W ) = 0.22264 

Gf = 1.16637- 10~ 5 GeV~ 2 aT 1 = 132.3407 . (3.1) 

We do not consider bottom and top quark effects. The remaining quarks are assumed to be 
massless and we work in the approximation where the CKM matrix is the identity matrix. 
We choose the invariant U W + W~^ mass as the central value for the factorization and 
renormalization scales: 



V-F = V-R = fio = Y (Pt! + Pi 2 + Pu! + Pu 2 + P-y) 2 - (3.2) 

We use the CTEQ6L1 parton distribution function at LO and the CTEQ6M set with 
a s (m z ) = 0.1176 at NLO [17]. 

Since the photon can be emitted either from the initial quark line or from each of the two 
final state charged leptons, we divide the phase space into three separate regions to consider 
all the possibilities and then sum the contributions to get the total result. The regions are 
generated as (approximately) on-shell W + — > e + z/ e 7 and W~ — > [i^Vp^ three-body decay 
and triple electroweak boson production, respectively. We impose a set of minimal cuts on 
the rapidity, y 7 (f>, and the transverse momenta , Pt jW , of the photon and the charged leptons 
which are designed to represent typical experimental requirements. Furthermore, leptons, 
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photons and jets must be well separated in the rapidity-azimuthal angle plane. Specifically, 
the cuts we impose are 

pr 7W >20GeV |y 7(£) |<2.5 R £l > 0.4 R je > 0.4 R h > 0.7 (3.3) 

where, in our simulations, a jet is defined as a colored parton of transverse momentum 
PTj > 20 GeV. For photon isolation, we implement the procedure introduced in [11]: if i 
is a parton with transverse energy Ej-. and a separation with a photon of transverse 
momentum p^ , then the event is accepted only if 

1 — COS 

E Ti 6(5 - R^) < p Tl (for all 5 < 5 ) (3.4) 

1 — cos do 

where 5q is a fixed separation that we set equal to 0.7. A quick look at Eq. (13.41) reveals that 
a sufficiently soft gluon can be arbitrarily close to the photon axis, while the energy of an 
exactly collinear parton must be vanishing in order to pass the isolation cut. Collinear-only 
events (leading to fragmentation contributions) are thus rejected, while soft emissions are 
retained, as desired. Additionally, for "ZZ7" production we impose that the invariant mass 
of any combination of two charged leptons, m«, be larger than 15 GeV in order to avoid 
virtual-photon singularities in 7* — > at a low m«. In Table [TJ we give results for the 
integrated cross sections for u W + W~ r j n and production at the LHC for the given 

cuts as well as for a harder cut on the photon transverse momentum of p^ > 30 GeV. In 
Table [21 the numbers for the Tevatron are presented for the cuts as in Eqs. (l3.3|l3.4p . but 
for less restrictive transverse momentum cuts, Pt iW > 10(10) GeV and Pt iW > 20(10) GeV. 
The center-of-mass energy is set to 14 TeV for the LHC and 1.96 TeV for Tevatron collisions, 
respectively. 



LHC 


LO [fb] 


NLO [fb] 


p T _ i(l) > 20(20) GeV 


1.695 


2.881 


p Ty(e) > 30(20) GeV 


0.9580 


1.738 


a("ZZ 7 " -> e + e-fi + fi--f) 
pr^ > 20(20) GeV 


0.07786 


0.1062 


a("ZZ 7 " -> e+e-/i + /i-7) 
p T ^ m > 30(20) GeV 


0.03969 


0.05577 



Table 1: Total cross sections at the LHC for pp — ► W + W~'y+X and pp -^■ZZ'y+X with 
leptonic decays, at LO and NLO, and for two sets of cuts. Relative statistical errors of the 
Monte Carlo are below 10~ 3 . 

In the following we show W + W~^ results for the process pp — > i/ e e + i)^~^ + X, i.e we do 
not consider interference effects due to identical leptons, which generally will be small. The 
cross sections have been multiplied by a combinatorial factor of 4, however, and correspond 
to the production of any combination of electrons or muons. 
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1 QlrOTTATl 

Ic VcLLIUIl 


t n [fKi 


nt n rfhi 


a( n WWi' -> e+i//i-z/ 7 ) 
■n^n ■> inn HpV 


0.9015 


1.271 


ff('¥^7" -> e+i//i-z/7) 
p T7W > 20(10) GeV 


0.3755 


0.5342 


a (»ZZi' -> e+e-/i+/i-7) 
pr, f<1 > 10(10) GeV 


0.05469 


0.07681 


a("ZZ 7 " -> e+e-/i + /i^7) 
Pr, f<1 > 20(10) GeV 


0.02440 


0.03441 



Table 2: Total cross sections at the Tevatron for pp — ► Vy + VF~7+V and — ►Z.Zt+JT with 
leptonic decays, at LO and NLO, and for two sets of cuts. Relative statistical errors of the 
Monte Carlo are below 10 -3 . 



In order to estimate the scale uncertainty of the total cross section, we study the variation 
of the numerical results in the interval 

Hf,Hr = £-Ho (0.1<f<10), (3.5) 

where /i is the li W + W~ n f" or u ZZ / y" invariant mass. The variation of the total u W + W~j" 
cross section at LO and NLO is shown on the left panel of Figure [2j It is evident that 
the LO scale variation greatly underestimates the size of NLO corrections, which give a PC- 
factor of 1.70 for n F = n R = /i . The NLO scale uncertainty is about 6% when varying the 
factorization and the renormalization scale \l — \ip — \ir up and down by a factor 2 around 
the reference scale /io — Qi slightly smaller than for the LO uncertainty, and is dominated 
by the dependence on /jr, which gives a negative slope with increasing energy, while it has 
a quite flat dependence on \i F . 

On the right panel of Figure (2J the scale dependence and relative size of the different 
NLO contributions are shown. One clearly sees that almost the entire scale variation of the 
integrated NLO cross section is accounted for by the real emission contributions, defined 
here as the real emission cross section minus the Catani-Seymour subtraction terms plus the 
finite collinear terms. The virtual contributions proportional to the Born matrix element 
(due to the M B term in Eq. fl2.ll) ) constitute the bulk of the NLO corrections, being roughly 
twice as large as the real emissions, while the finite virtual remainders due to boxes and 
pentagons only represent 3% and less than 1% of the total result, respectively, and their 
scale dependence is basically flat. 

Similar considerations apply for "ZZ7" production where a combinatorial factor of 2 has 
been included for pp — > ZZj + X — > ifi^i^i^l + X. In Figure [3], one observes that the 
overall scale variation decreases when varying the factorization and the renormalization scale 
A* = = ^0 U P ari d down by a factor 2 around the reference scale /i = Q from roughly 
± 7% at LO to ±1% at NLO, and we have a K-factor of 1.36 for Hf = = Ho- The 
main difference with respect to the u W + W~ , ~y" case is the numerical contribution of the real 



7 





solid: H F = H R = t, |i 


15 






dashed: n F = S M-o . m= 






14 


dotted: n F = Ho. ^R= 5^0 




— Total NLO 


12 




10 


Virtual-born 


10 

8 


NLO 
LO 


D 

5 


Real 








Virtual-box 


6 


i 





Virtual-Pentagons 

i 



0.1 1 10 0.1 1 10 



Figure 2: Left: Scale dependence of the total LHC cross section for pp — ► W + W~'y + X — ► 
£ + £~'-f + jr3 T + X at LO and NLO within the cuts of Eqs. I\3.3i3.4\) . The factorization and 
renormalization scales are together or independently varied in the range from 0.1 -p, to 10 -hq . 
Right: Same as in the left panel but for the different NLO contributions at fi F = /ir = £/ir> 



corrections, which amounts to 20 % of the total contribution or less. 

The NLO contributions cannot be described by a simple rescaling of the LO results since 
the corrections show a strong dependence on the phase space region under investigation. As 
practical examples, we plot two differential distributions at LO and NLO together with the 
associated K-factor, defined as 

da NLO /dx 

= ^JclV ■ (3 - 6) 

The observables, x, considered are the transverse-momentum of the photon and the separa- 
tion between the softest (lowest-pr) lepton and the photon in W + W~ r )' production at the 
LHC. In Figure H] we show the transverse-momentum distribution of the on-shell photon: 
the corresponding K-factor increases from 1.5 at pT-y = 20 GeV to 2.1 at pt-j = 140 GeV. In 
Figure we consider the separation between the photon and the softest lepton and a similar 
variation of the K-factor in the range 1.5 — 1.9 is observed. 



4 Conclusions 



We have calculated the NLO QCD corrections to the processes pp,pp — > W + W~ , y+X 
and pp,pp — ► ZZ'-f+X with full leptonic decays of the W and Z bosons. These processes 
can be relevant both for New Physics searches (as a background) and for the measurement 
of quartic gauge couplings (as a signal) at the LHC. 

Our numerical results show sizeable cross section increases with respect to the LO cal- 
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Figure 3: Left: Scale dependence of the total LHC cross section for pp — > ZZ'-f + X — > 
ifi^i^i^l + X at LO and NLO within the cuts of Eqs. ( j3.3f3.4j) . Tie factorization and 
renormalization scales are together or independently varied in the range from 0.1 -u. to 10- Hq. 
Right: Same as in the left panel but for the different NLO contributions at p, F = fi R = £/x . 




20 40 60 80 100 120 140 20 40 60 80 100 120 140 

Figure 4: Left: Transverse-momentum distribution of the photon in W + W~^ production at 
the LHC. LO and NLO results are shown for \ip = u.r = /x and the cuts of Eq. (j3.3j) . Right: 
K-factor for the transverse- momentum distribution of the photon as defined in Eq A3.6\) . 
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Figure 5: Left: Separation between the photon and the softest lepton in W + W~^ production 
at the LHC, at LO and NLO with \ip — [ir = yUo and the cuts of Eq. I\3.3\) . Right: K-factor 
for the separation between the photon and the softest lepton as defined in Eq. A3.6\) . 

culation, ranging from 40% to as much as 100% in certain regions of phase space. W + W~ r y 
and ZZ^ production at the LHC provide additional examples of cross sections whose the- 
oretical errors at LO are substantially underestimated by considering scale variations only: 
the LO factorization scale variation is much smaller than the NLO correction. Remaining 
NLO scale variations are at the 10% level (e.g. ±6% for the integrated W + W~ r y production 
cross section at the LHC when varying /ir = u. F = u. by a factor of 2 around the reference 
scale fi = Q = m WWl ). 

Given the size of the higher-order corrections and, in particular, their strong dependence 
on the observable and on different phase space regions under investigation, a fully-exclusive 
NLO parton Monte Carlo for W + W~ r ) and ZZ^ production is required to match the ex- 
pected precision of the LHC measurements. We plan to incorporate this and other processes 
with a final state photon into the VBFNL0 package in the near future. 
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